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We discuss various analytical approximation methods for the evolution of the 
density fluctuation in the Universe. From primordial density fluctuation, the large- 



scale structure is formed via its own self-gravitational instability. For this dynamical 
evolution, several approaches have been considered. In Newtonian cosmology, in 
which matter motion is described by Newtonian dynamics with cosmic expansion, 

Lagrangian description for the cosmic fluid is known to work rather well for quasi- 

>' 

■ nonlinear clustering regime. In this paper, we briefly review Lagrangian perturbation 

(N ' 

. theory in Newtonian cosmology. 
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£j : I. INTRODUCTION 

c3 



There are various structures in the Universe. For example, galaxies, group of galaxies, 
cluster of galaxies, voids, large-scale structure, and so on. How are such structures formed? 
As the most plausible scenario, we have considered that a primordial density fluctuation 
spontaneously produces by itself a gravitational instability. The fluctuation could have 
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34, 351. When we do 



been generated by several processes in the early Universe 
not consider the superhorizon scale or extremely dense region like a supermassive black 
hole, the motion of the cosmological fluid can be described by Newtonian cosmology. In 
Newtonian cosmology, the motion of matter is described by Newtonian dynamics, and the 
cosmic expansion is given by the Friedmann equation. 

The Lagrangian description for the cosmological fluid can be usefully applied to the 
structure formation scenario. This description provides a relatively accurate model even in 



a quasi-linear regime. Zel'dovich |63( proposed a linear Lagrangian approximation for dust 
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fluid. This approximation is called the Zel'dovich approximation (ZA) P. llll 0. 0. I 
After that, hig her-order approximation for the Lagrangian description was proposed 
fxoL Il2 . [T3L Il4 0. 0]. Then it was shown that the Lagrangian approximation describes 
the evolution of density fluctuation better than the Eulerian approximation 000. 
Especially for the planar model, ZA gives an exact solution 0,0]. 

Although the Lagrangian approximation gives an accurate description until a quasi-linear 
regime develops, it cannot describe the model after the formation of caustics. For example, 
in ZA, even after the formation of caustics, the fluid elements keep moving in the direction 
set up by the initial condition. Therefore, the nonlinear structure that forms diffuses. 

In order to proceed with a hydrodynamical description in which caustics do not form, 
a qualitative pressure gradient pj] and thermal velocity scatter 0, in a collisionless 



was 



matter have been discussed. Additionally the "adhesion approximation" (AA) 
proposed based on the nonlinear diffusion equation (Burgers' equation). In AA, an artificial 
viscosity term is added to ZA. Because of the viscosity term, we can avoid caustics formation. 
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Using AA, the problem of structure formation has been discussed [31, 
density divergence does not occur in AA, and the density distribution close to the N-body 
simulation can be produced. However, the origin of the viscosity has not yet been clarified. 
Another modified model for the Lagrangian approximation is known as the "Truncated (or 
Optimized) Zel'dovich approximation" (TZA) 0. I^. The evolution equation in TZA is 
the same as that in ZA. Instead of the evolution, the initial power spectrum is changed 
in TZA. Because the caustics are formed by small-scale fluctuation, TZA introduces some 
cutoff in the small scale in the initial spectrum. As a result, TZA avoids the formation of 
the caustic. However, the physical meaning of the cutoff also has not yet been clarified. 

For another purpose, improvement of approximation, several modified models have been 
proposed. For example, Bagla and Padmanabhan p] developed a new scheme that extended 
the frozen flow approximation 39]. Then they showed better agreement with the N-body 
simulations than ZA. Instead of ordinary perturbative expansion, the Pade approximation 
was introduced to Lagrangian approximation 0, . When we consider the spherical void 
evolution, if we increase the order of Lagrangian approx imation to try to improve, contrary 
to expectation, the description becomes worse [4J, |54(. On the other hand, the Pade ap- 
proximation improves the description of the void evolution. Another approach is generalized 
ZA. Although ZA solutions are independent of position, the solutions of generalized models 
depend on position. Although the computation becomes complicated, we can obtain exact 
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solutions for planar, cylindrical, and spherical symmetric models |^, 13, |37 

In ZA and its extended models, pressure was ignored. Recently, Lagrangian approxi- 
mation in which the effect of pressure was taken into consideration have been analyzed. 

discussed the effect of velocity dispersion using the collision- 
They argued that models of a large-scale structure should be 



Buchert and Dominguez 
less Boltzmann equation 
constructed for a flow describing the average motion of a multi-stream system. Then they 
showed that when the velocity dispersion is regarded as small and isotropic it produces ef- 
fective "pressure" or viscosity terms. Furthermore, they derived the relation between mass 
density p and pressure P, i.e., an "equation of state." Buchert et al. 3] showed how the 
viscosity term or the effective pressure of a fluid is generated, assuming that the peculiar 
acceleration is parallel to the peculiar velocity. Dommguez clarified that a hydro- 

dynamic formulation is obtained via a spatial coarse-graining in a many-body gravitating 
system, and that the viscosity term in AA can be derived by the expansion of coarse-grained 
equations. 



With respect to the relation between the viscosity term and effective pressure, and the 
extension of the Lagrangian description to various matter, the Lagrangian perturbation the- 
ory of pressure has been considered. Adler and Buchert |2| have formulated the Lagrangian 
perturbation theory for a barotropic fluid. Morita and Tatekawa 43J and Tatekawa et al. 
57| solved the Lagrangian perturbative equations for a polytropic fluid up to the second 
order. Hereafter, we call this model the "pressure model." The behavior of the pressure 
model strongly depends on the polytropic exponent 7 0, I&Z, Q According to a recent 



study 



the effect of pressure cannot realize the artificial viscosity in AA completely. 



Matarrese and Mohayaee proposed Lagrangian perturbative equation for two component 
(baryon and cold dark matter) fluids . They consider lowest-order Lagrangian perturba- 
tion for the two components and discussed the evolution of the density fluctuation. 

We especially mention only various Lagrangian perturbation models. For advanced and 
related topics, we recommend to refer other reviews and books (0 0, 0, for 
example). 

This paper is organized as follows. In Sec. |Hj we present basic equations for the cos- 
mo logical fluid in Newtonian cosmology. First, we show Eulerian description for the basic 
equations in Sec. Ill Al Then we introduce Lagrangian description in Sec. Ill Bl In Sec. lIIII we 
show perturbative solutions for dust fluid up to a third-order approximation. First, we con- 
sider first-order longitudinal solutions in Sec. IIII Al The first-order approximation is called 
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"Zel'dovich approximation" (ZA). After that, a higher-order approximation is proposed. In 
Sec. IIIIBI we mention second- and third-order longitudinal solutions. In Sec. IIII CI we ex- 
plain the transverse mode. When we analyze vorticity for cosmological fluid, we need to 
consider this mode. Although these solutions seem useful, they do not have physical mean- 
ing after formation of a caustic. To avoid formation of a caustic, several modified models 
have been proposed. In Sec. IIII Dl and IIII El we show modified models for Lagrangian per- 
turbation. In Sec. IIII Dt we explain some modification to avoid the formation of caustics. 
In Sec. IIII El we show alternative models to obtain high accuracy. In Sec. IIVI we explain 
the pressure model. In Sec. we summarize the success and failure in Lagrangian approx- 
imation and state applications. We notice the spherical void evolution with dust fluid in 
Appendix El In Appendix^ we provide a the table explaining the physical quantities used 
in this paper. 



II. BASIC EQUATIONS 

A. Basic equations for the cosmological fluid 

In this section, we present basic equations for the cosmological fluid in Newtonian cosmol- 
ogy; i.e., the motion of fluid is described by Newtonian dynamics, and the background cosmic 
expansion is given by the Friedmann equation. In this situation, the motion of perfect fluid 



is described by a continuous equation, Euler's equation, and Poisson's equation 
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|f) +V r -(pu) = 0, (1) 



r 



{ du\ 1 
- + ( tt .V) rU =--V r P + g, (2) 

g = r -V r $ ,V • g = AnGp , (3) 

To introduce cosmic expansion, we adopt the coordinate transformation from physical co- 
ordinates to comoving coordinates. 

T 

x = —r^r, a(t) : scale factor , (4) 
a(t) w w 

where r and x are physical coordinates and comoving coordinates, respectively. The scale 
factor satisfies the Friedmann equations: 

'a\ 2 9 8nG JC A 
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a AttG A , . 

a = -— pt+ 3' (6) 

H^-, (7) 

with a curvature constant /C and a cosmological constant A. H = a/a and are Hubble 
parameter and background density, respectively. Here we define density parameter Q: 

Under the transformation (j3J), the velocity and partial differential operator are changed as 

u = r = ax + v(x, t), (v = ax) , (9) 

V x = aV r , (10) 
( df(x = r/a,t) \ fdf\ a 

Here we define the density fluctuation 5 as follows: 

p = p b (t){l + 6(x,t)} , p b (xa- 3 . (12) 
In the comoving coordinates, the basic equations for cosmological fluid are described as 

! + V W l + *)} = 0, (13) 

dv 1 . „ , a 1 1 „ „ . „ . 

_ + -(«. V x )v + -v = -g V X P , (14) 

at a a a ap 

V x xg = 0, (15) 
V x g = -4nGp b a6. (16) 



Therefore, we obtain from Eqs. (fT3j) and (|T5|) 

ut a a 



v - x h^7 + + - ( u - v -H= ' ( 17 ) 



and from Eqs. (fUjl and (fTTIj) . 

V - • (iST + + " ( v ■ V -) v l = ~^Gp h 5 - V x ■ (—V x p) . (18) 
\ot a a J \pa J 

B. The Lagrangian description for the cosmological fluid 

In the Lagrangian hydrodynamics, the coordinates x of the fluid elements are represented 
in terms of Lagrangian coordinates q as 

x = q + s(q,t), (19) 
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where q is defined as initial values of x, and s denotes the Lagrangian displacement vector 
due to the presence of inhomogeneities. The peculiar velocity is given by 



v = as . 



(20) 



Then we introduce the Lagrangian time derivative: 



d 9 1 „ 

at at a 

Using the Lagrangian derivative, we rewrite Eqs. (JT7j) and (JTSJ). The nonlinear term of the 
peculiar velocity disappears. 

„ / dv a \ . . 

v * x (* + ^) = ' (22) 

V, ■ + %) = -4*G Pb 6 - V, • Q^pj . (23) 

The exact form of the energy density in the Lagrangian space is obtained from Eq. (j!3|) 

as 

P = Pb J- 1 , (24) 

where J = det (dxi/dqj) = det(5ij+dsi/dqj) is the Jacobian of the coordinate transformation 
from x to q. J is described by expansion of the derivative of Lagrangian perturbation as 
follows: 

J = 1 + V... + i((V f ..)'-gQ + d*(g). (25) 

Next, we consider the description of the vorticity in Lagrangian space. We define the vorticity 

as 

lj = -V x xv. (26) 

(X 

From Eq. (122(1 . we obtain the vorticity equation. 

duj h u? . . _ 

+ + - (V, •«) = («• V x x) . (27) 
dr a a 

In the Lagrangian description, Eq. (|27j) is solved exactly. 



, , (go(g) • V 9 )a;(g,t) 

W = ^7 ' (28) 



where cjo(q) is the primordial vorticity. 

From Eqs. (J2Uj) and (|23jh the peculiar gravitational field is written as 

* = "( i+2 |*-?fw j " lv - J )- < 29 » 
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where an overdot (') denotes d/dt. Hence, from Eqs. f)15jl and (fTfij). we obtain the following 
equations for s: 

V* x (s + 2-s) = 0, (30) 



a 

V x ■ (s + 2-J - r l V x J^j = -4 7 rGp b (J- 1 - 1) . (31) 

If we find solutions of Eqs. ()30|) and ()31|) for s, the dynamics of the system considered is 
completely determined. Since these equations are highly nonlinear and hard to solve exactly, 
we will advance a perturbative approach. Remark that, in solving the equations for s in 
the Lagrangian coordinates q, the operator Va; will be transformed into V q by the following 
rule: 

d dxj d d 8sa d 

= I = | L [32) 

dqi dqi dxj dxi dqi dxj 

We decompose s into the longitudinal and the transverse modes as s = V g S + S T with 
V q -S T = 0. 

III. DUST MODEL 
A. Zel'dovich approximation 



Zel'dovich derived a first-order solution of the longitudinal mode for dust fluid [63[. The 
evolution equation for first-order solution is written as 



£(i) + 2 -S^ - 47rGp b S w = . (33) 
a 

Then, the first-order solutions are written as follows: 

SW(q,t) = D + {t)^\q) + D4t)^\q) . (34) 

For the E-dS Universe model (Q = 1), D± is written as 

D+(t) = t 2/3 oc a(t) , (35) 
D_(t) = r 1 oc a(t)- 3 / 2 . (36) 

For a non-flat model, following Doroshkevich et al. , we replace the standard cosmological 
time by a new time r defined by 

dr oc a~ 2 dt . (37) 
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Using new time, we obtain an analytic form for D±. For the open model (Q < 1), 

D + (r) = l + 3 (r 2 -l)(l + Ilog(^)) , (38) 
D_(r) = r(r 2 -l). (39) 

For the closed model (Q > 1), we obtain the analytic form by the transformation r — > ir in 
Eq. (jSi). 

For A-flat model, in which the cosmological constant exists and the Universe is flat, D± 
is written as 

poo AQ h /5 2\ 

D +M = */» = ^Me's)' (40) 

D-(t) = h, (41) 
H = %, (42) 

where Bx/ h 2 is incomplete Beta function: 

B z (fx,u)= f Z -pr\l- p y- l 6.p. (43) 
Jo 

This first-order approximation is called the Zel'dovich approximation (ZA). Especially when 
we consider the plane-symmetric case, ZA gives exact solutions P. I^. 

Using Eqs. pijl and (jHlJ), we can describe density fluctuation. Here we define the defor- 
mation tensor X a p\ 

X a/3 = ^ = 6 aP + ^. (44) 
dq/3 dq 



The eigenvalues of the deformation tensor (}44j) are written by 

Wi = -D(t)\°(q) , (45) 

where A° is the eigenvalue of dtp^/dqp. Using these eigenvalues, the density fluctuation is 
described by 

(1 + J)" 1 = Eti(l + «;*)■ (46) 

When Wi becomes —1, the caustic ( "Zel'dovich's pancake") is formed and the density fluc- 
tuation diverges. After that, the perturbative solution does not have physical meaning. 
In order to avoid the formation of caustics, several modified models have been proposed. 
We will discuss these models in Sec. IIII Dl By generalizing the eigenvalues of the deforma- 
tion tensor (Eq. 1)450) . more exact approximations can be proposed. These models will be 
expressed in Sec. IIII El 
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B. Higher-order approximation 



ZA solutions are known as perturbative solutions, which describe the structure well in 
the quasi-nonlinear regime. To improve approximation, higher-order perturbative solutions 
of Lagrangian displacement were derived. Irrotational second-order solutions (PZA) were 

n n 

derived by Bouchet et al. M and Buchert and Ehlers [13J], and third-order solutions (PPZA) 
were obtained by Buchert jlJ|. Bouchet et al. ^(|, and Catelan jl7|. 

Hereafter we consider only growing-mode D + in first-order Lagrangian perturbation. 
First, we show the evolution equation for second- and third-order perturbative equations 
for the longitudinal mode. Here we decompose the Lagrangian perturbation to time and a 
spacial component. 



S®(t,q) = E(t)^(q), 
S®(t,q) = F(t)^(q), 



(47) 
(48) 



where the superscript means n-th order Lagrangian perturbation. The second-order 
perturbative equation is written as 



E + 2-E - ^Gp b E) 

a J ' 

-2nG Pb D 2 MST-i>Ml 



(49) 



The third-order perturbative equation becomes little complicated. 



F + 2-F - AnGp b F i/j 



,(3) 



= -8nG Pb [D{E - D 2 ) (^VS " 1>%W$) + D 3 det(^) 
The spacial component of the second order are written as follows: 

^ = (^S) 2 -^g. 



The time component of second-order E obeys the evolution equation. 



E + 2-E - 4nGp b E = -2irGp b D z 
a 



(50) 



(51) 



(52) 



For the third order, we divide to two components of which one is derived from combining 
^M 1 ) and i})( 2 ' and the other is derived from (-?// 1 ') 3 . 



S {3) = F a (t)^\q) + F b (t)^ 3 \q), 
ifr® u = det , 



(53) 
(54) 
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*ft-*sws-*sw5. (55) 



a 



F a + 2-F a - AnGp b F a = -MGp b D A , (56) 



a 

2\ 



F b + 2-F b - 4irGp b F b = -8nGp b D(E - D 2 ) . (57) 

(X 

For the E-dS Universe model, E and F and described with analytic form. 

E{t) = - fr* x D{t) 2 , (58) 
F a (t) = -lt*<xD{t) 3 , (59) 
F b (t) = fe<xD{ty. (60) 

For the non-flat Universe model, E can be described with analytic form. For the open 
Universe model, 



*<r) = -i-!<r" 



r f t — 1\ 1 / r 2 - 1 /r - 1 

1+ 2 log 7TT + 2 T + — log 0m 



(61) 



where r is defined by Eq. ()37|) . For the closed model (fi > 1), we obtain analytic form by 
the transformation r — > ir as in Eq. (|38j) . 

For the A-flat model, unfortunately, we need to solve evolution equations with the nu- 
merical method [10]. 

3(h 2 - l)E + 2hE = 2E - 2D 2 , (62) 

3(h 2 - l)F a + 2hF a = 2F a - 2D 3 , (63) 

3(h 2 - l)F b + 2hF b = 2F b + 2D 3 (l - . (64) 

C. The transverse mode 

For the transverse mode, until third-order solutions have been obtained (H Q, 

Q- The 

evolution equation for the first-order perturbation is given as 

S T + 2*S T = 0. (65) 
a 

For the E-dS model, the first-order solution is written as 

S T ( q ,t) = Sj(q) + t-^Snq). (66) 



Therefore the transverse mode does not have a growing solution in the first-order pertur- 
bation. Furthermore, if we consider only the longitudinal mode for the first order, the 
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second-order transverse mode does not appear. However, the third-order transverse mode 
is generated by the triplet of the first-order longitudinal mode, and it grows. Although 
the third-order solution has the growing mode, this solution does not show that a vorticity 
grows. 



D. Modified models I — avoidance of formation of a caustic 



Cosmological N-body simulations show that pancakes, skeletons, and clumps remain dur- 
ing evolution. However, when we continue applying the solutions of ZA, PZA, or PPZA after 
the appearance of caustics, the nonlinear structure diffuses and breaks. 

Adhesion approximation (AA) 27] was proposed from a consideration based on Burgers' 
equation. This model is derived by the addition of an artificial viscous term to ZA. AA with 
small viscosity deals with the skeleton of the structure, which at an arbitrary time is found 
directly without a long numerical calculation. 

We briefly describe the adhesion model. In ZA, the equation for "peculiar velocity" in 
the E-dS model is written as follows: 

du 



da 



+ (u ■ V x )u 







u 



dx 

da 



x 

a 



(67) 
(68) 



where a(oc t 2 / 3 ) is scale factor. To go beyond ZA, we add the artificial viscosity term to the 
right side of the equation. 



du 

da 



+ (u ■ V x )u = uVlu. 



(69) 



We consider the case when the viscosity coefficient v — > +0 [y ^ 0). In this case, the 
viscosity term especially affects the high-density region. Within the limits of a small v, the 
analytic solution of Eq. (JBTIj) is given by 



u (x 



= E 



x - q a 



Ja exp ( ) / Y^Ja eX P ( 



2v 



(70) 



where q a is the Lagrangian points that minimize the action 



I a = I(x, a; q a ) = S (q a ) + 



x 



2a 



mm. 



Ja 

So 



det ( 5ij + 



d 2 S 
dqidqj 



-1/2 



q=Qo 



S(q,t c 



(71) 
(72) 
(73) 
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considered as a function of q for fixed x [3J|. In AA, because of the viscosity term, the caustic 
does not appear and a stable nonlinear structure can exist. However, when the evolution 
is advanced too much, the Universe becomes covered with a high-density structure called 
"skeleton." 

As another modification th eory , truncated (or Optimized) Zel'dovich approximation 



(TZA) has been used well 



421 ]. During evolution, the small scale structure contracts 



and forms caustics. Therefore, if we introduce some cutoff in the small scale, we will be able 
to avoid the nation of ca^icHQ. 

In TZA, for the avoidance of caustics, we introduce a cutoff in the initial density spectrum. 
Various cutoff methods were considered, and the Gaussian cutoff was the most suitable from 
the comparison with N-body simulation. The Gaussian cutoff is introduced to the initial 
density spectrum as follows: 

V(k, t in ) — ► V(k, t- m ) exp (~k 2 /k NL ) , (74) 

where k^L is the "nonlinear wavenumber", defined by 

_ r^NL 

l = D + (t) 2 V{k,t- m )dk, (75) 



'ko 

where D(t)+ is the growing factor in ZA (Eq. (JSJ)). Because the "nonlinear wavenumber" 
depends on the time, the wavenumber becomes small during evolution. Therefore, when the 
evolution is advanced too much, small structure will vanished. 



E. Modified models II — high-accuracy models 

Bagla and Padmanabhan developed a new approximation scheme that improved ZA. 
They considered extension the frozen flow approximation 3t|. Then they showed better 
agreement with the N-body simulations than ZA. 

n 

Yoshisato et al. [62] introduced the Pade approximation for the evolution of the density 
fluctuation. The Pade approximation seems to be a generalization of the Taylor expansion. 
For a given function f(x), the Pade approximation is written as the ratio of two polynomials: 

m * 1 +h ak f k , (76) 

i + Efc=i b kx k 

where a k and b k are constant coefficients. Assume we already know the coefficient q (0 < 
/ < M + N) of the Taylor expansion around x — 0. 

M+N 

f(x)= E c lX l + o{x M+N+1 ). (77) 
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Comparing the coefficents a^, b k , and c&, we determine a k and b^. 

a = c , (78) 
= h mCk~m (k = l,---,N), (79) 

m=l 

N 

b m c N _ m+k = -c N+k (k = 1, • • • , N) . (80) 

m=l 

The advantage of the Pade approximation is that even if we consider a same-order expansion, 
the Pade approximation describes original function rather better than the Taylor expansion. 

First, Yoshisato et al. [62] showed that the application of the Pade approximation for 
Eulerian description improves the past Eulerian models. Then Matsubara et al. proposed 
the application of the Pade approximation for Lagrangian description and improved PZA 
and PPZA. 

Another approach is known as "local" approximations. In ZA, the growing factor D is 
independent of the position. On the other hand, in local approximations, the factor depends 
on position. 

First, we derive an evolution equation for the eigenvalue of the deformation tensor 
(Eq. (jlll ). Here we use time variable r (Eq. Ij37|) ). Differentiating of Eq. (f2T)|) with re- 
spect to Xj and summing up the trace part, we obtain 

X ^^f = -^Ga 4 p b (J- 1 - 1) . (81) 



If is diagonal, 



Eq. (jHTj) is written as 



X tj = (1 + w^Sij , (82) 



i (1 + Wi) 



AttGcl 4 



Pb 



(83) 



(l + wi) (1+^2X1+^3) 
For a spherical case we have wi = u?2 = ^3; for a cylindrical case w± = W2 and W3 = 0; and 
for a planar case w<i = W3 = 0. ZA is a linearized form for Wi of Eq. (|83p. 

Reiseneger and Miralda-Escude proposed changing the eigenvalue (Eq. (|%5)l) as 

w % — ~D(t, A°)A°(q) . (84) 

This is then substituted in Eq. (|83jl . They named this model the modified Zel'dovich ap- 
proximation (MZA) or extension of ZA (EZA). The MZA is exact for planar, spherical, and 
cylindrical symmetric cases. However, they also reported that the MZA may not work for 
underdense regions. 
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Audit and Alimi .4] introduced another ansatz. Eq. (JHHj) may be written in this form: 

3 



2^ (1 + w j +Wk + WjW k )Wi - AnGa p b I H 1 — | ir, 



i=i 



0. 



(85) 



They split this equation into three equations for each wf. 



[1 + Wj + Wk + WjWk)Wi = AirGa p& I 1 + 



Wj + Wk WjWk 



Wi . 



(86) 



This approximation is called the Deformation Tensor Approximation (DTA) Although 
this is also exact for planar, spherical, and cylindrical symmetric cases, the splitting of 
Eq. ()85)1 is not unique, and we could add more local terms in Eq. (|86J) . 

Betancort-Rijo and Lopez- Corredoida |7] proposed another generalization. In terms of 
the linear solution Aj = D(t)\®(q) in ZA, they generalized Eq. (|4*3j) to 



Wi = -ri(\i,\j,\k)\°i(q) 



(87) 



where rj(Aj, A,-, A^) is the power series of A. 



r j (A i> A i> A fc ) = l+ E Cl nhn {\ 3 + \k)\\ 3 -\k? n X 

Lm,n=0 



in 
i i 



(88) 



where Cf m n are the coefficients of the p-th order terms (p = l + 2n+m). The ZA corresponds 
to Ti = 1. They named this generalized model the "Complete Zel'dovich approximation" 
(CZA). They calculated explicitly the coefficients Cf mn up to the fourth order of A in the 
E-dS Universe model 3). This model is also exact for planar, spherical, and cylindrical 
symmetric cases. Furthermore, the CZA describes the evolution of ellipsoid dust better 
than the MZA, the DTA, and ZA [37]. However the CZA does not apply for perturbations 
with negative values of A° (?[ 37 1. 



IV. PRESSURE MODEL 



Although AA seems a good model for avoiding the formation of caustics, the origin of the 
modification (or artificial viscosity) is not clarified. Buchert and Dominguez [15J argued that 
the effect of velocity dispersion becomes important beyond the caustics. They showed that 
when the velocity dispersion is still small and can be considered isotropic, it gives effective 
"pressure" or viscosity terms. Buchert et al. [16] showed how the viscosity term is generated 
by the effective pressure of a fluid under the assumption that the peculiar acceleration is 
parallel to the peculiar velocity. 
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Adler and Buchert j^J have formulated the Lagrangian perturbation theory for a 
barotropic fluid. Morita and Tatekawa 4^ and Tatekawa et al. [57] solved the Lagrangian 
perturbation equations for a polytropic fluid in the Friedmann Universe. Hereafter, we call 
this model the "pressure model" . 

If the pressure is given by a barotropic equation of state, the pressure affects only the 
longitudinal mode in the first-order. Therefore, the transverse solution in the pressure model 
becomes the same as that in the dust model. The evolution equations are written as 



5 T + 2-5 T = 0, (89) 



a 



S + 2-S- A7rGp b S - 4 ^ 
a a 2 dp 



V 2 S = 0. (90) 



p=pb 



When we assume the polytropic equation of state P = np 1 , the first-order equation for 
the longitudinal mode is written as follows. 

S + 2-S - AnG Pb S - ^l^V 2 S = . (91) 
a a 2 

The solution manifestly depends on the scale of the fluctuation. We adopt Fourier transfor- 
mation in Lagrangian space for Eq. (j91jl . 

S + 2-S -A7TGp b S + K1Pb \K\ 2 S = , (92) 
a a 2 

where K is a Lagrangian wavenumber. In the E-dS Universe model, we can write the 



first-order solution relatively simply |43j. For 7 7^ 4/3, 

S(K, a) oc a- 1 '* J ±5/{8 -^ {{^^^ ^ )/2 ) > (93) 
where J v denotes the Bessel function of order u, and for 7 = 4/3, 

S(K, a) oc a -i/4±V25A6-c 2 |K|W ; (94) 

where C% = 47rGpb(Qin) a in/3 and C 2 = ^7Pb(0'm) 7_1 a in^ 7_1 ^- a 'm is scale factor when an 
initial condition is given. 

For the other FRW universe model, if 7 takes a special value, we obtain an analytic 
solution with a hypergeometric function j^: For open or closed model, 7 = 1,4/3. For the 
A-flat Universe model, 7 = 1/3,4/3. 

In this model, the behavior of the solutions strongly depends on the relation between the 
scale of fluctuation and the Jeans scale. Here we define the Jeans wavenumber as 

K >=[dPlMri ■ (95) 
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S(K,a) oca 1/4 J± 5 /(8-6 7 ) 



(96) 



43 



The Jeans wavenumber, which gives a criterion for whether a density perturbation with a 
wavenumber will grow or decay with oscillation, depends on time in general. If the polytropic 
index 7 is smaller than 4/3, all modes become decaying modes and the fluctuation will 
disappear. On the other hand, if 7 > 4/3, all density perturbations will grow to collapse. 
In the case where 7 = 4/3, the growing and decaying modes coexist at all times. 
We rewrite the first-order solution Eq. ()93|) with the Jeans wavenumber: 

In the pressure model, the second-order perturbative solutions have been derived 
Because of mode coupling, the calculation of the second-order perturbative solutions are 
complicated. Only for the E-dS Universe model and the case where 7 = 4/3 have the 
solutions of the longitudinal and the transverse mode been derived Then, in the E- 

dS Universe model, if 7 is larger than 4/3, it was shown that the effect of the second-order 
pressure becomes weak gradually during evolution [5j| . In this paper, the authors considered 
the expansion of the Bessel function, 

Jv{Z) -\2) h^ + n + l)- (97) 

In the E-dS Universe model, z oc a^ 4 " 37 ^ 2 (Eq. (J93|) ) . If 7 is larger than 4/3, the index of 
argument of the Bessel function becomes negative. Therefore, precision does not become 
much worse even if we drop the higher terms of the expansion. Using finite expansion, they 
analyzed the second-order perturbation for the case of 7 = 5/3. 

These models are considered for a single component fluid. Of course, there are several 
components in the Universe, for example, baryonic matter and dark matter. Matarrese and 
Mohayaee 



38j proposed the Lagrangian perturbation model for two component fluids. One 
is dark matter which affects only gravity. Another is baryonic matter which affects gravity 
for all fluid and pressure for baryonic matter. 

When we consider more than one species of matter, the evolution equation of the density 
fluctuation is modified 48j. We assume a polytropic equation of state for baryonic matter. 

P oc p 7 . (98) 

Taking the divergence of the baryon Euler equation, we obtain the evolution equation for 
the baryon density fluctuation with the Lagrangian perturbation. 

3 

'^B^ — ^T. TTTTTV^ 1 ^°b) 



(7-l« 



(99) 
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where B and DM mean baryon and dark matter, respectively, kj is the Jeans wavenumber 
of the baryonic matter. (') denotes the derivative by a scale factor (' = d/da). Here we 
assume that the amount of the baryonic matter is much less than that of the dark matter. 
Therefore, we ignore the gravitational force of the baryonic matter. Substituting Eq. (|25j) 
for (jnnj), we obtain 

' + b- ' '» - ~ir^ - w^m v2 ( W ■ (100) 

Next, we consider only first-order terms for the Lagrangian displacement s. 



V, " s" B + ^V, ■ s' B - ^V q ■ s DM = ■ s B . (101) 



We note that the same Eulerian position x(t) is generally reached by the two components 
from different Lagrangian positions q B and qoM- Here we ignore this difference and set 
simply qs = qDM, because we consider the lowest order approximation. 

We assume that the perturbation has only a longitudinal mode (s = V q S). Under 
this assumption, we obtain the evolution equation of the Lagrangian perturbation for the 
baryonic matter. 

3 3 1 3 

S * + Ya S ' B ~ 2o^¥^ Sb = Ya S ° M ■ (1 ° 2) 

This equation can be solved in Lagrangian Fourier space. 



V. SUMMARY 



We briefly review Lagrangian perturbation theory in Newtonian cosmology. From a com- 
parison to Eulerian perturbation theory, Lagrangian perturbation theory describes quasi- 
nonlinear evolution of density fluctuation rather well. We notice that when we analyze phys- 
ical phenomena or statistical quantities with Lagrangian perturbation, because the physical 
coordinates are Eulerian ones, we must transform the coordinates from Lagrangian ones to 
Eulerian ones. Although transformation from Eulerian ones to Lagrangian ones is given by 
Eq. (|19jl . the inverse transformation seems a little difficult. 

Although Lagrangian perturbation theory seems to work well for describing quasi- 
nonlinear evolution of the density fluctuation, it is not always useful. For example, we 
consider the evolution of the spherical void (4^, |54j . In general, when we consider higher- 
order perturbation, the approximation is improved. However, for void evolution, when we 
advance for a long time, higher-order terms show bad behavior. The even-order (2nd, 4th,...) 
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perturbative terms promote to development of the reverse, i.e., the void disappears and the 
density fluctuation converges. On the other hand, the odd-order (3rd, 5th,...) perturbative 
terms promote evolution excessively. We provie details in Appendix 1X1 



In Sec. IIIID[ we noticed modified models for avoidance of formation of a caustic. However, 
the origin of those modifications is not clarified. Buchert et al. [16] discussed the relation 
between the viscosity term in the Adhesion approximation and the effective pressure of 
a fluid. Can the pressure become the origin of the modification? According to a recent 
study the pressure term can realize that the density distribution looks like that of AA 
until reaching quasi-nonlinear regime. However, because of Jeans instability, the behavior 
of the density fluctuation in the pressure model seems different to that in AA. 

On the other hand, the pressure especially affects small-scale structure. Therefore, we can 
expect that the pressure term realize the Gaussian cutoff in TZA. In fact, the correspondence 
between "nonlinear wavenumber" k^i in TZA (Eq. (175(1 ) and the Jeans wavenumber Kj in 
the pressure model (Eq. (j95(l ) has been discussed [5i|. As a result, the correspondence has 
not been explained sufficiently The reason is as follows: First, in TZA, k^L affects only 
the initial spectrum. On the other hand, Kj affects the evolution of fluctuation. Second, 
although fc^vL obviously depends on the initial spectrum, we did not clarify the dependence 
on the initial condition of Kj. We think that a consideration of the physical process, which 
was not considered here, or the analysis of the N-body simulation, is necessary for a decision 
of Kj, i.e. k. 

Therefore, the effect of pressure cannot explain modifications in the Lagrangian approxi- 
mation. If we consider the origin of these modification with the effect of velocity dispersion, 
we may have to consider anisotropic velocity dispersion 36^ . 

The Lagrangian perturbation theory seems rather useful until quasi-nonlinear regime de- 
velops. However, the density fluctuation becomes strongly nonlinear in the present Universe. 
Can we apply the theory to the problem of structure formation? Because a huge simulation 
has executed ((3])' someone may thinks that the perturbation theory is no longer useful. 
We think that analytic methods are important for examining the physical process of struc- 
ture formation. For example, in the discussion about statistical quantities, the analytical 



method shows remarkable results 
One possibility is analysis of t 
projects have been progressing |l 



20, 



45 



49|,Eil|- 



re past structure. Recently, several galaxy redshift survey 
2fiJ, |2£j. In these projects, many galaxies within a region 



(z < 0.3 for 2dF, z < 0.5 for SDSS, and 0.7 < z < 1.4 for DEEP2) have been observed. 
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For the next generation of spectroscopic survey, for example, Kilo-Aperture Optical Spec- 
trograph (KAOS) has been proposed 30(. The two primary scientific purposes of KAOS are 
(i) the determination of the equation of state of dark energy and (ii) the study of the origin 
of galaxies. As one of the other science applications for KAOS, the growth of structure is 
considered. In this project, many high-z galaxies will be observed. Therefore, this project 
would not only probe topology of the matter distribution but also the dynamical history of 
structure formation. Because the density contrast will still be small in the high-z region, 
we expect that we will be able to discuss the characters of the density fluctuation using the 
Lagrangian perturbation theory well. 

Recently, various dark matter models have been proposed j^J. Some of them affect not 
only the gravity but also a special interaction. We also show that the linear approximation of 
the pressure model seems rather good until a quasi-linear regime develops. If the interactions 
of the dark matter are affected by effective pressure, the linear approximation can be applied 
for the analysis of the quasi-nonlinear evolution of the density fluctuation. 
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APPENDIX A: SPHERICAL VOID EVOLUTION WITH DUST FLUID 



We describe the development of the spherical void with Lagrangian perturbation. Here 
we consider "top-hat" spherical void, i.e., a constant density spherical void. In the E-dS 
Universe model, the equation of motion of a spherical shell is written as 



d / 2 dx 
dt\~dt 



2a 2 x 



X J 



where x is a comoving radial coordinate and xq 
5 = a for a — > 0, Eq. (jAljl can be integrated. 

2 



(Al) 



x(to) |4jJ. Under the initial condition 



'cLffY 
da 



1 3 



(A2) 



where R{9) = a(t)x/xo is physical particle trajectory. The exact solution for the expansion 
of a top-hat void (Eq. (|A2JI ) can be parameterized as follows: 

3 



R{9) 
a(9) 



10 

3 
5 



(cosh 9 - 1] 
- (sinh 9 — 



l 2/3 



(A3) 
(A4) 



From these equations, we can obtain density fluctuation. 

3 



5(x) 



X ) 

9(e-sinh^ 2 



- 1 



2(cosh#- l) 3 

The perturbative solution for the void in the E-dS Universe model is given by 



(A5) 



R(t) = Ro 



1 - J2(-l) k C k a k 



k=l 



(A6) 



where C k is given by Lagrangian perturbative equations. Munshi, Sahni, and Starobin- 

n 

sky 44] derived up to the third-order perturbative solution (Ci,C 2 ,C 3 ). In addition to 
these, Sahni and Shandarin ^ obtained C4 and C5. They concluded that ZA remains the 
best approximation to apply to the late-time evolution of voids. Especially, if we stop to 
improve until even order (2nd, 4th, 6th, ■ ••), the perturbative solution describes the con- 
traction of void. From the viewpoint of the convergence of the series, the conclusion seems 
reasonable. Table H] shows the coefficients of Lagrangian perturbation theory for top-hat 
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TABLE I: The ratio of coefficients of Lagrangian perturbation theory for top-hat spherical voids 
(Cfc_x/Cfc). The higher order the perturbation we consider, the worse the approximation becomes in 
late-time evolution. In other words, when the order of the perturbation increases, the convergence 
of the series becomes bad. 



k 2 a 3 a A b 5 b 6 

Ck-x/Ck 7 3.52174 2.80517 2.49236 2.31645 



k 



10 



11 



Ck-i/Ck 



2.20357 



2.12495 



2.06704 



2.02260 



1.98743 



"Munshi, Sahni, and SlcCTbbinsky 
b Safmi and Shandarin [54J 



spherical voids. Here we derived up to the eleventh order. The convergence radius of the 
series a c is defined by the ratio of the coefficients: 

a c = lim — — . (A7) 

When the scale factor a becomes larger than a c , the series (Eq. (|A6p ) does not converge. 

We estimate the convergence radius from finite series. Table H] shows the ratio Ck-i/Ck- 
When k is increased, C^-i/Ck decreases. Therefore, the higher order the perturbation we 
consider, the worse the approximation becomes in late-time evolution. If we consider only 
normal Lagrangian perturbation, the problem cannot be solved. The Pade approximation 
[40! lo^ shows that the approximation for the late-time evolution of voids is improved, and 
therefore some applications are expected. 



APPENDIX B: QUANTITIES USED IN THIS PAPER 

In this paper, we use the following variables for physical quantities, coordinates, and so 

on. 



TABLE II: Coordinate and velocity variables. 



r 


physical coordinates 


X 


comoving Eulerian coordinates 


Q 


Lagrangian coordinates 


u 


velocity in physical coordinates 


V 


peculiar velocity 







TABLE III: Time variables. 



t 


standard cosmological time 


r 


new time variable (Eq. (|87|)) 


a 


scale factor 


z 


red shift 



TABLE IV: Cosmological parameters. 



H = a/a 


Hubble parameter 


K 


curvature constant 


A 


cosmological constant 


n 


density parameter 



TABLE V: Physical quantities of matter. 


P 


matter density 


Pb 


background (averaged) matter density 


P 


pressure of matter 


G 


gravitational constant 


$ 


gravitational potential 


9 


gravitational force in comoving coordinates 


6 


density fluctuation (Eq. (|12j0 


s 


Lagrangian displacement (perturbation) vector 




vorticity (Eq. (J27J)) 



TABLE VI: Lagrangian perturbation. 


S 


Laerranenan nerturbation notcntial florientiiHiTial mode) 


T 

S 


Lagrangian perturbation (transverse mode) 


J 


Jacobian of the coordinate transformation from x to q 




spacial part of Lagrangian perturbation potential 




eTowiner factor in ZA (Fin i 

w i w v v ill ^ iUi\_/ uv/i ill i—i -i \ J — J v_| • tM^^^^j! 1 


D- 


decaying factor in ZA (Eq. (1341)) 


E 


prnwinpr fartnr in PZA fFVi (I4QD^ 


F 


growing factor in PPZA (Eq. (|5U|)) 




deformation tensor (Eq. (|44jl) 


A? 


eigenvalue of dib^/dqp 


Wi 


eigenvalue of X a p (Eq. (H5J)) 



TABLE VII: Some quantities appearing in the pressure model. 



7 

K 

K 

Kj 
kj 



polytropic index in equation of state 
proportional constant in equation of state 
Lagrangian wavenumber vector 

Jeans wavenumber in Lagrangian coordinates (Eq. ()95jl ) 
Jeans wavenumber in Eulerian coordinates 



